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Executive  Summary 


This  report  is  part  of  an  ongoing  effort  to  develop  computer  models  of  human  response  to 
combined  injuries  anticipated  after  a  nuclear  detonation.  The  following  study  describes  im¬ 
provements  and  extensions  of  a  mathematical  model  of  the  small  intestine  epithelium.  The 
model  is  designed  to  simulate  small  intestine  crypt  and  villus  cellular  populations  following 
exposure  to  acute  radiation  and  thermal  injuries.  This  information  can  be  used  to  quantify 
damage  of  the  epithelial  lining  of  the  small  intestine. 

Major  results  of  this  modeling  effort  include: 

•  A  murine  model  has  been  optimized  to  a  large  set  of  radiation  response  data. 

•  A  human  model  has  been  optimized  to  radiation  response  data. 

•  Burn  effects  of  the  murine  model  have  been  optimized  against  thermal  injury  data. 

•  Burn  response  of  the  human  model  has  been  extrapolated  from  the  murine  model. 

•  The  murine  model  has  been  validated  with  radiation  and  burn  combined  injury  data. 
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1  Introduction 

Applied  Research  Associates,  Inc.  (ARA)  has  been  tasked  by  the  Defense  Threat  Reduction 
Agency  (DTRA)  to  support  its  mission  to  safeguard  the  United  States  against  weapons  of 
mass  destruction  (WMDs).  ARA  supports  this  effort  by  developing  state-of-the-art  mathe¬ 
matical  models  that  predict  medical  and  performance  effects  of  radiation  and  burn  injuries, 
thereby  enhancing  our  understanding  of  the  potential  impact  of  a  nuclear  detonation. 

1.1  Modeling  Gastrointestinal  Cell  Kinetics 

Radiation  and  burn  injuries  affect  the  cellular  composition  of  the  epithelial  lining  of  the  small 
intestine.  A  reduced  epithelial  lining  results  in  nutrient  absorption  insufficiencies  and  per¬ 
meability  changes.  Understanding  the  cellular  kinetics  of  small  intestine  epithelial  cells  will 
aid  in  predicting  and  analyzing  casualties.  To  this  end,  we  developed  mathematical  models 
of  the  cellular  kinetics  of  the  small  intestine  in  response  to  radiation,  burn  and  combined 
injuries.  Murine  models  have  been  constructed  to  assist  with  human  model  development. 

1.2  Model  Implementation 

The  models  developed  in  this  report  are  being  implemented  in  HENRE  (Health  Effects  from 
Nuclear  and  Radiological  Environments),  a  user-friendly  software  platform  with  which  one 
can  (1)  define  an  insult,  (2)  run  models  to  predict  health  effects  of  that  insult,  and  (3) 
analyze  model  output. 

1.3  Existing  Mathematical  Models 

There  are  a  number  of  mathematical  models  of  the  small  intestine  that  quantify  the  effects 
of  radiation.  Of  these  models,  two  stand  out  in  particular.  Russian  scientist  Olga  Smirnova 
has  developed  a  mathematical  model  that  simulates  three  cellular  compartments  of  the  small 
intestine  epithelium:  proliferating  crypt  cells,  maturing  crypt  cells  and  villus  cells  (Smirnova, 
2010).  The  model  is  structured  as  a  system  of  ordinary  differential  equations  (ODEs)  with 
nonlinear  proliferation  and  feedback  effects.  Radiation  is  assumed  to  damage  a  portion  of 
the  proliferative  cells  in  this  model  using  the  one-target-one-hit  theory  of  cell  damage  (Lea, 
1955) 

Another  powerful  model  of  small  intestine  epithelial  response  is  built  into  the  Radiation- 
Induced  Performance  Decrement  (RIPD)  model  (Anno  et  ah,  1989,  Anno  et  ah,  1991). 
RIPD,  a  physiologically-based  mathematical  model  of  the  effects  of  ionizing  radiation,  calcu¬ 
lates  time-dependent  symptoms  and  predicts  performance  following  exposure  to  irradiation. 
The  gut  injury  model  (GIM),  which  is  built  into  RIPD,  measures  the  effects  of  ionizing  ra¬ 
diation  on  the  small  intestine  epithelial  lining.  This  model  is  nested  into  three  components: 
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the  lethal  and  potentially  lethal  (LPL)  model  (Curtis,  1986)  calculates  probabilities  of  sur¬ 
vival  for  irradiated  cells;  the  proliferation  and  intracellular  repair  (PAIR)  model  determines 
proliferation  rates  and  mitotic  delay  from  the  LPL  outputs;  and  a  compartmental  model 
determines  the  number  of  clonogenic  cells,  transit  cells  and  villus  cells  from  the  LPL. 

Compared  to  Smirnova’s  model,  the  GIM’s  quantification  of  DNA  damage  and  DNA 
repair  is  more  detailed  than  the  cellular  damage  component  of  Smirnova’s  model.  This 
detail  is  particularly  useful  for  simulating  protracted  exposure,  as  the  time  components  of 
damage  and  repair  are  crucial  for  capturing  dose  response  for  various  dose  rates.  We  are 
currently  concerned  with  the  effects  of  acute  radiation  where  exposure  time  is  generally  low, 
and  dose  rate  is  generally  high,  and  thus  Smirnova’s  model  of  cellular  damage  is  sufficient. 
The  GIM  is  currently  built  into  HENRE  and  we  can  thus  leverage  its  capabilities  for  future 
projects  that  may  involve  simulating  intestinal  response  to  protracted  radiation  exposures. 

Smirnova’s  model  is  more  detailed  than  the  GIM  at  the  cellular  compartment  level,  as 
it  includes  more  interactions  and  feedback  effects.  This  granularity  improves  the  recovery 
kinetics  of  the  model,  which  are  important  for  accurate  assessment  of  an  exposed  individual. 
Smirnova’s  model  also  leverages  more  biological  data  to  support  values  chosen  for  model 
parameters.  We  have  made  improvements  to  Smirnova’s  model  for  our  study  which  include 
fitting  the  model  to  a  large  set  of  experimental  data,  and  adding  burn  response,  which  is 
absent  in  Smirnova’s  original  model. 

There  are  additional  models  of  the  small  intestine  epithelium,  each  with  downfalls  and 
limitations.  As  discussed  in  Smirnova,  2010,  these  include  a  model  which  is  unable  to 
accurately  capture  cellular  damage  from  radiation  (Meineke  et  ah,  2001),  a  model  which 
requires  a  heavy  amount  of  data  to  accurately  capture  radiation  response  (Tyazelova,  1988), 
and  a  model  which  is  too  complex  to  achieve  reasonable  parameter  identification  (Gozenbuk 
and  Keirim-Markus,  1988). 
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2  Methods 

The  following  section  provides  details  of  onr  parameterization  approach,  including  cost  func¬ 
tion  definitions  and  the  optimization  algorithm  used  for  parameter  fitting.  We  used  different 
approaches  for  parameterizing  the  human  and  murine  models  because  there  is  a  large  amount 
of  rodent  data  and  very  little  human  data.  The  purpose  of  constructing  the  murine  model 
is  to  assist  the  development  of  the  human  model  through  species-to-species  parameter  ex¬ 
trapolations. 

2.1  Murine  Model  Parameterization 

Here,  we  describe  our  approach  for  parameterizing  the  murine  small  intestine  model.  We 
initially  found  biological  measurements  of  mice  including  cell  counts,  migration  rates,  cell 
cycle  times,  etc.  Accounting  for  variability  in  the  data,  we  set  ranges  for  parameters  directly 
dependent  on  these  measurements.  This  is  discussed  in  more  detail  in  Section  3.5. 

After  exhausting  available  biological  measurements  to  inform  parameter  bounds,  we  fit 
model  parameters,  excluding  those  involved  in  burn  response,  to  25  experimental  time  series 
datasets.  This  data  includes  21  radiation  response  datasets  (Section  3.5.2)  and  four  datasets 
of  response  to  insults  other  than  radiation  or  burn  (Section  B.l).  All  parameters,  exclud¬ 
ing  those  involved  in  burn  response,  were  fit  to  this  data  using  an  optimization  algorithm 
described  in  Section  2.3.  During  optimization,  parameter  bounds  for  biologically  supported 
parameters  were  enforced,  and  large  ranges  were  allowed  for  parameters  lacking  prior  knowl¬ 
edge.  The  following  cost  function,  Cm,  was  used  to  determine  the  difference  between  murine 
model  simulations  with  a  parameter  set  p,  and  the  25  time  series  datasets: 

25 

Cm(p)  =  '£E,(p).  (2.1) 

i=l 

The  error  from  simulating  each  time  series,  Eu  is  calculated  as  an  average  of  absolute 
residuals  between  the  time-dependent  experimental  data  (  (ti,di),  ...,  ( tn,dn )  )  and  the 
trajectory  of  the  corresponding  simulated  variable  (v(t)): 

1  n 

E{p)  =  -J2\di~v(ti)\-  (2-2) 

n7E[ 

We  use  the  sum  of  the  absolute  residuals,  as  opposed  to  the  more  conventional  sum  of  squared 
errors,  to  measure  error  in  this  study  for  multiple  reasons.  For  one,  defining  the  cost  function 
with  the  sum  of  squared  errors  assumes  that  the  errors  are  normally  distributed,  which  is 
not  necessarily  true  here.  Also,  the  data  is  diverse  and  is  expected  to  contain  outliers  due  to 
the  diversity  of  experimental  procedures.  For  instance,  the  data  used  for  the  murine  model 
are  from  various  sections  of  the  small  intestine  and  multiple  species  (rats  and  mice).  By 
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using  the  absolute  residuals,  we  place  less  emphasis  on  outliers  of  the  data  than  we  would 
with  squared  errors. 

Burn  response  was  implemented  in  the  murine  model  by  adding  functions  constructed 
from  empirical  evidence.  More  detailed  biological  data  is  required  to  model  this  response 
mechanistically.  The  parameters  added  to  the  model  from  these  functions  were  fit  to  exper¬ 
imental  time  series  data  using  the  same  optimization  algorithm  used  to  fit  the  first  set  of 
parameters.  Much  less  data  was  available  to  fit  the  burn  parameters  (2  time  series  datasets, 
as  opposed  to  25). 

2.2  Human  Model  Parameterization 

Similar  to  our  approach  with  the  murine  model,  we  first  used  experimental  measurements 
to  justify  boundaries  for  as  many  parameters  in  the  human  model  as  possible  (see  Section 
3.5  for  details). 

We  defined  the  cost  function  differently  for  the  human  parameterization  due  to  the  fact 
that  we  only  found  one  time  series  dataset  which  we  could  simulate  with  the  model.  The 
cost  function  for  the  human  model  Ch,  dependent  on  the  parameter  set,  p,  compares  the 
normalized  epithelial  data  from  Trier  and  Browning,  1966,  where  13  epithelial  measurements 
are  taken  over  the  course  of  radiation  treatment  (  (ti,di),  ...,  ( tm,dm ),  m=13  )  to  z(t), 
the  simulated  normalized  villus  cell  counts.  This  cost  function  places  extra  emphasis  on 
simulating  the  nadir  of  the  data  set,  which  occurs  at  (t$,d$): 


Ch(p)  =  e(p )  +  En{p),  (2.3) 

where  errors  from  simulating  the  average  data  point  (E)  and  the  nadir  (En)  are  defined 
below: 


i  m 

E(p)  =  ~  “  di\,  (2-4) 

En(p)  =  \z(t8)  -  d8\.  (2.5) 

We  place  particular  emphasis  on  the  nadir  of  the  human  data,  which  we  calculate  in  both 
of  the  expressions  above,  and  thus  has  weight  We  believe  the  nadir  of  the  cell  counts 
to  be  an  important  biomarker  for  the  health  of  an  exposed  individual. 

2.3  Optimization  for  Parameter  Fitting 

When  searching  for  an  optimal  parameter  set  to  best  fit  data,  we  solve  the  model  using  the 
“ode”  function  of  the  deSolve  package  (Soetaert  et  ah,  2015)  in  the  statistical  computing 
environment  R  (R  Core  Team,  2013).  We  define  the  cost  function  using  the  modCost 
tool  in  the  Flexible  Modelling  Environment  (FME)  library  (Meysman  et  ah,  2008),  and  we 
perform  optimization  using  the  DEoptim  tool  in  the  Global  Optimization  by  Differential 
Evolution  library  (Ardia  et  ah,  2015).  DEoptim  uses  the  standard  search  algorithm  known 
as  differential  evolution  (DE),  to  search  parameter  space  for  a  global  optimal  parameter  set. 
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3  Modeling  the  Small  Intestine 

This  section  provides  biological  and  mathematical  details  of  the  cellular  kinetic  model  of  the 
small  intestine.  Following  a  description  of  the  model,  we  present  the  data  simulated  by  the 
model  and  provide  simulation  results  of  radiation  and  burn  experiments. 


3.1  Biological  Description 


The  small  intestine  connects  the  stomach  and  the  large  intestine  in  the  gastrointestinal  tract. 
The  small  intestine  is  divided  into  three  segments:  the  duodenum,  the  smallest  segment 
which  connects  to  the  stomach;  the  jejunum,  the  largest  segment;  and  the  ileum,  the  final 
section  connecting  the  jejunum  to  the  large  intestine. 

The  epithelial  lining  of  the  small  intes¬ 
tine  is  composed  of  projections  called  villi. 

Multiple  compartments  located  around  the 
base  of  each  villus,  called  crypts,  supply  the 
villi  with  epithelial  cells.  In  humans,  about 
five  crypts  supply  each  villus  (Ferguson  et 
ah,  1977),  while  there  are  about  10  crypts 
per  villus  in  mice  (Wright  and  Irwin,  1982). 

Throughout  the  following  report,  we  refer 
to  the  collection  of  a  villus  and  its  associ¬ 
ated  crypts  as  a  column.  Clonogenic  cells 
located  near  the  base  of  the  crypts  provide 
the  column  with  new  cells,  which  proliferate 
as  they  migrate  up  the  crypts  (Quastler  and 
Sherman,  1959).  Proliferating  cells  begin  to 
mature  as  they  reach  their  respective  villus. 

Once  on  the  villus,  the  cells  (now  referred  to 
as  villus  cells)  travel  to  the  villus  tip,  where 
they  are  shed  into  the  lumen  of  the  intestine 
(Figure  3.1). 

Radiation  and  burn  injuries  have  an  im¬ 
mediate  effect  on  crypt  and  villus  epithe¬ 
lial  cells  in  the  small  intestine.  Proliferating 
crypt  cells  respond  to  radiation  by  immedi¬ 
ately  initiating  a  mitotic  delay,  and  many  proliferating  cells  undergo  apoptosis  (Lesher,  1967, 
Bond  et  ah,  1965).  Burn  has  an  immediate  effect  on  villus  cells,  promoting  early  cell  death 
before  migration  is  complete  (Carter  et  ah,  2014,  Wolf  et  ah,  1999).  This  response  is  fast 
and  strong  but  is  not  observed  after  two  days  in  mice  (Wolf  et  ah,  1999).  Crypt  cells  exhibit 
an  additional  response  to  burn  which  dampens  proliferation.  This  secondary  response  has 


Figure  3.1:  Diagram  of  crypt  and  villus  struc¬ 
ture  and  epithelial  cell  dynamics. 
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been  reported  in  mice  no  sooner  than  one  day  after  burn,  but  persists  for  at  least  a  week 
post-insult  (Jeschke  et  al.,  2007). 


3.2  Mathematical  Model 

The  following  section  describes  the  mathematical  model  of  the  small  intestine  we  have  de¬ 
veloped  for  this  study.  Figure  3.2  provides  a  diagram  of  the  model  including  the  maturation 
and  migration  of  the  small  intestine  epithelial  cells,  feedback  elements  of  the  system,  and 
the  effects  of  radiation  and  burn.  Negative  feedback  through  the  proliferation  compartment 
from  crypt  and  villus  cells  is  modeled  based  on  the  chalone  feedback  reported  in  numerous 
cell  lines  by  Bullough,  1971.  In  Smirnova,  2010,  this  feedback  is  assumed  to  depend  on  a 
mediator,  promoted  by  crypt  and  villus  cells,  which  inhibits  crypt  cell  proliferation.  Here, 
we  modify  this  assumption,  proposing  that  a  generic  mediator  stimulates  proliferation,  but 
the  level  of  this  mediator  is  down- regulated  by  the  crypt  and  villus  cells.  We  have  made  this 
modification  in  anticipation  of  simulating  drug  applications  which  generally  promote,  rather 
than  inhibit,  cellular  processes.  Changing  this  assumption  does  not  have  a  significant  effect 
on  the  structure  of  the  mathematical  model.  In  addition  to  negative  feedback  through  prolif¬ 
eration,  there  is  also  evidence  that  maturation  rates  through  crypt  and  villus  compartments 
are  directly  related  to  the  rate  of  proliferation  (Smirnova,  2010).  We  will  discuss  responses 
to  radiation  and  burn  in  Sections  3.3  and  3.4,  respectively. 


Cell 

transitions 

* 

Transitions 
affected  by 
radiation 

I 

Transitions 

affected 
by  burn 

-  > 

Feedback 

Effects 
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Healthy 

cells 
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by  radiation 
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Mediator 
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Cell  Loss 

c 

Cell 
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Figure  3.2:  A  diagram  of  the  small  intestine  cell  kinetic  model  with  burn  and  radiation 
responses. 
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The  complete  mathematical  model  of  the  small  intestine  and  a  comprehensive  list  of 
variables,  parameters,  and  parameter  values  are  provided  in  Appendix  A.  The  model  is  a 
set  of  ordinary  differential  equations,  following  the  structure  in  Smirnova,  2010,  tracking 
proliferating  crypt  cells  (x),  maturing  crypt  cells  (y)  and  villus  cells  (z),  each  representing 
the  number  of  cells  per  villus  column.  In  the  model,  the  generic  mediator  /  is  assumed  to 
upregulate  cellular  proliferation  (Equation  3.7).  At  the  same  time,  each  cell  type  is  assumed 
to  downregulate  the  mediator  (Equation  3.4).  The  transit  rates  F  for  maturing  crypt  cells, 
and  E  for  villus  cells  are  assumed  to  be  linearly  dependent  on  the  proliferation  rate  5,  as 
seen  in  Equations  3.5  and  3.6.  This  is  based  on  experimental  evidence  that  suggests  both 
independent  and  dependent  transit  rates  for  both  cell  types  (Smirnova,  2010). 

Equations  of  the  model,  excluding  radiation  and  burn  effects,  are  provided  below: 


x  =  Bx  —  'yx 

(3.1) 

V  =  ix-  Fy 

(3.2) 

z  =  Fy  —  Ez 

(3.3) 

i  =  G  -  H(x  +  9y  +  -&z)I 

(3.4) 

F  —  5(1  +  LB) 

(3.5) 

E  =  0(1  +  MB) 

(3.6) 

B=  a 

1  +  K/I 

(3.7) 

In  the  above  equations,  the  following  are  constant:  7,  G,  H,  6 ,  $,  5,  L,  -0,  M,  a  and  K. 
Assuming  mediator  interactions  occur  on  a  fast  time  scale,  we  solve  for  /  using  a  quasi- 
steady-state  approximation.  Setting  /  =  0  in  Equation  3.4  gives  us  the  following: 


H(x  +  6y  +  ftz) ' 

Substituting  Equation  3.8  into  Equation  3.7, 


(3.8) 


where 


B 


a 

1  +  P(x  +  0y  +  •dz) 


P 


HK 

~G~ 


(3.9) 

(3.10) 


Compared  to  the  model  structure  of  Smirnova,  2010,  it  is  clear  that  changing  the  mediator 
effect  from  inhibitory  to  stimulatory  only  affects  the  definition  of  the  constant  /3. 

At  equilibrium,  x  =  y  =  z  =  0,  x  =  x,  y  =  y  and  z  —  z.  Assuming  x  >  0,  y  >  0,  and 
z  >  0,  the  following  equations  represent  the  relative  equilibrium  compartment  sizes  (using 


Equations  3.1  -  3.3): 


a/7  —  1 

/3(1  +  5(1+L7)  +  V>(l+Af7) ) 

g  =  7 

"0(1  +  M7) 

We  use  normalized  numbers  of  crypt  and  villus  cells  to  simulate  normalized  experimental 
data.  For  this  reason,  we  normalize  the  variables  of  the  mathematical  model,  which  we 
express  as  x ,  y  and  z.  A  non-dimensionalization  of  the  model  variables  can  be  found  below: 


(3.11) 

(3.12) 

(3.13) 


7x(.B  —  1) 

(3.14) 

5  ((1  +  l)x  -  (1  +  lB)y) 

(3.15) 

/  ^  _U  jji  _  _  \ 

0(^1  +  /(1  +  lB)y  (l  +  mB)zj, 

(3.16) 

where  l  =  7 L,  m  =  7 M  and  B  =  B/ 7. 

When  simulating  experimental  data,  we  use  z  to  represent  the  normalized  number  of 
villus  cells  per  column.  We  simulate  normalized  crypt  cells  per  column  using  the  crypt  cell 


mass  (CCM  =  —  — ) . 
v  x  +  yJ 


This  can  be  expressed  in  terms  of  the  normalized  variables: 


CCM 


5(1  +  l)x  +  7  y 
5(1  +  /)  +  7 


(3.17) 


Similarly,  we  can  express  the  normalized  total  number  of  epithelial  cells,  or  epithelial  mass 
x  y  ~\~  z 

(EM  =  - - - - ).  in  terms  of  normalized  variables: 

v  x  +  y  +  z' 


_  05(1  +  0(1  +  m)x  +  70(1  +  m)y  +  57(1  +  l)z 

05(1  +  /)(1  +  m)  +  70(1  +  m)  +  57(1  +  l) 


(3.18) 


3.3  Radiation  Injury 

Following  exposure  to  radiation,  proliferating  crypt  cells  become  damaged  and  die  by  apop¬ 
tosis  (Potten  and  Grant,  1998).  Following  the  structure  of  the  model  in  Smirnova,  2010,  we 
account  for  this  damage  by  dividing  the  proliferation  compartment  into  damaged  ( xd )  and 
undamaged  ( xud )  cellular  sub-compartments  (x  =  xud  +  xd).  We  assume  an  immediate  cel¬ 
lular  exchange  from  undamaged  to  damaged  compartments,  simulating  dose  response  with 
a  multitarget  single- hit  model  (Joiner,  2009).  This  theory  proposes  that  one  hit  of  radiation 
in  n  different  sensitive  DNA  regions  (targets)  is  required  to  kill  a  radiosensitive  cell.  In 
response  to  a  dose  of  D  Gy  at  t  —  0,  we  set  initial  values  of  xud  and  xd  to  the  following: 


9 


D/D°\n 


(3.19) 

(3.20) 


a0(O)  =  x{\  —  e 
xud(0)  =x(l  -xd(0)) 

where  the  parameters  D0  and  n  determine  the  degree  of  the  dose  response.  We  assume  the 
damaged  cells  die  at  an  average  rate  of  is,  but  these  cells  can  feed  back  through  a  mediator 
to  inhibit  proliferation.  0  represents  the  strength  of  this  feedback,  relative  to  the  strength 
of  feedback  by  undamaged  cells.  Updating  the  rate  of  proliferation  to  reflect  this: 


B 


a 

1  +  /3(xud  +  (j)xd  +  9y  +  ftz)  ’ 


and  normalizing  (xud  =  xud/x  and  xd  =  xd/x ),  we  find  the  following: 


B  = 


ah 


1  +  b(xud  +  (j)xd  + 


7 


+  T&Z 


7 


5(1  +  1))  0(1  +  m) 


) 


where 


b 


a/'y  —  1 
$7  $7 

5(1  +  Z)  0(1  +  m) 


(3.21) 


(3.22) 


(3.23) 


3.4  Thermal  Injury 

There  is  very  little  data  on  human  small  intestine  epithelial  cell  response  to  burn,  but  these 
effects  have  been  documented  in  rodents.  Multiple  studies  of  mice  and  rats  have  reported  that 
villus  cells  are  quickly  killed  through  apoptosis  immediately  following  burn  (Carter  et  ah, 
2014,  Chen  et  ah,  2013,  Jeschke  et  ah,  2007,  Ramzy  et  ah,  2000,  Wolf  et  ah,  1999).  In  mice, 
this  response  is  strong  and  rapid,  generally  lasting  less  than  two  days  (Wolf  et  ah,  1999). 
In  addition  to  villus  killing,  there  is  a  delayed  downregulation  of  crypt  cell  proliferation. 
In  Jeschke  et  ah,  2007,  a  60%  total  body  surface  area  (TBSA)  burn  in  rats  resulted  in 
significant  proliferation  suppression  that  lasted  for  more  than  7  days,  while  proliferation  was 
observed  to  be  normal  one  day  after  burn.  We  incorporate  empirical  functions  in  our  model 
to  account  for  these  effects  and  fit  the  parameters  of  these  functions  to  experimental  data. 

In  our  model,  we  simulate  response  to  a  burn  of  S%  TBSA  by  modifying  the  rates  of 
villus  cell  loss  and  cellular  proliferation  with  functions  W0(S,t)  and  W\(S,t),  respectively: 


xud  =(B~  Wi(S,  t))  xud  -  ~ixud  (3.24) 

z  =  Fy-(E  +  Wo(S,t))z  (3.25) 


where 
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(3.27) 


Wi(S,t) 


{s 

bi  —  —  -  exp  (  -  ai£(100  -  5))  £  >  tb  +  Ab 

0  t  <  4  +  Ab 


In  Equations  3.24  -  3.27,  0  <  S'  <  100,  4  is  the  time  of  burn  onset  and  A*,  is  the  delay- 
after  burn  before  proliferation  suppression  begins,  a 0  and  cq  determine  the  duration  of 
the  burn  effects,  b0  and  b\  determine  the  maximum  effects  of  burn,  and  k  is  an  activation 
threshold  for  the  burn  effect.  Assuming  an  equivalent  proliferation  suppression  delay  in  mice 
and  rats,  we  set  At  =  1  day  in  the  murine  model,  reflecting  the  results  reported  in  Jeschke 
et  ah,  2007. 


Wq  and  W\  have  similar  structures,  where  bt 


S 

S  +  k 


i 


0, 1  represents  the  initial  response 


to  burn,  and  exponential  functions  are  used  to  represent  the  relaxation  of  the  burn  responses. 
W\  employs  an  additional  dependence  on  burn  size  in  the  time  component  of  the  relaxation 
function  to  account  for  the  fact  that  proliferation  suppression  lasts  much  longer  for  larger 
burns  (Jeschke  et  ah,  2007)  than  smaller  burns  (Baker  and  Valeriote,  1968,  Carter  et  ah, 
2013,  Wolf  et  ah,  1999). 


3.5  Experimental  Data 

The  following  section  presents  the  data  used  in  this  study  to  parameterize  and  validate  the 
human  and  murine  small  intestine  models.  As  discussed  in  Section  2,  our  general  approach  for 
model  parameterization  is  the  following:  utilize  biological  measurements  of  epithelial  cells 
to  inform  as  many  parameters  as  possible,  fit  remaining  unknown  parameters  (excluding 
those  associated  to  burn  response)  to  acute  radiation  data  and  fit  the  burn  parameters  to 
acute  burn  data.  We  have  no  human  burn  data,  so  we  chose  to  extrapolate  the  human  burn 
parameters  from  the  murine  model. 

3.5.1  Biological  Data  Used  to  Inform  Parameters 

We  have  been  able  to  leverage  biological  measurements  of  epithelial  cell  compartments  and 
kinetics  to  inform  a  subset  of  parameters  for  the  murine  and  human  models.  Table  3.1 
lists  biological  measurements  of  epithelial  cells  in  the  murine  small  intestine,  and  Table  3.2 
provides  the  relationships  between  these  measurements  and  the  model  parameters.  Using 
the  information  from  Tables  3.1  and  3.2,  we  placed  bounds  on  murine  model  parameters  7, 
a,  5,  and  i/j.  We  used  the  same  approach  for  the  human  model,  where  Tables  3.3  and  3.4  list 
human  biological  measurements  and  their  relationships  to  model  parameters,  respectively. 


Tabic  3.1:  References  and  bounds  for  cellular  measurements  in  the  murine  small  intestine. 


Value  measured 

Intestinal 

Section 

Reference 

Lower 

Upper 

(Eq)  Transit  rate 
through  villus 
compartment  (d-1) 

Jejunum 

Wright  and  Irwin,  1982 

0.506* 

0.83* 

(CP)  Crypt 

population 

(cells/crypt) 

Jejunum 

Wright  and  Irwin,  1982 

367.7* 

525* 

{GF)  Growth 
fraction 

Jejunum 

Al-Dewachi  et  al.,  1975 

0.613° 

0.65° 

{CCPD)  Cells 
produced  in  crypt 
(cells/d) 

Jejunum 

Wright  and  Irwin,  1982 

254* 

394* 

{MIC)  Crypt 
mitotic  index  {%) 

Ileum 

Matsuzawa  and  Wil¬ 
son,  1965 

2.88** 

2.88** 

- 

Smirnova,  2010 

5.3°° 

5.3°° 

(Tm)  Length  of  M 
phase  (h) 

D  uodenoj  ej  unum 

Wright  et  al.,  1973 

1.09+ 

1.31+ 

(7s)  Length  of  S 
phase  (h) 

Jejunum 

Schultze  et  al.,  1972 

O 

8° 

*Lower  or  upper  bound  for  95%  Cl. 

**Only  one  measurement  provided. 

°Two  values  were  calculated  using  different  measuring  techniques. 

00  Original  reference  unknown. 

1  Human  measurements,  only  two  provided. 

Table  3.2:  Murine  parameter  associations  with  experimental  measurements. 


Parameter 

Biological  Meaning 

Formula^ 

7 

Steady-state  rate  of  proliferative  compartment 

MIC 

repopulation 

Tm  *  GF 

a 

Maximum  rate  of  proliferative  compartment 
repopulation 

1/  (7m  +  7s) 

s 

Minimum  transit  rate  through  maturation 

CCPD 

compartment 

CP{1-GF){1  +  1) 

Minimum  transit  rate  through  villus  compart¬ 
ment 

E0/{  1  +  m) 

§  Experimental  measurements:  crypt  mitosis  index  (MIq),  growth  fraction  ( GF ), 
crypt  cells  produced  per  day  ( CCPD ),  crypt  population  {CP),  duration  of  M 
phase  (7m)  and  duration  of  S  phase  (7s). 
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Table  3.3:  References  and  bounds  for  cellular  measurements  in  the  human  small  intestine. 


Value  measured 

Intestinal 

Section 

Reference 

Lower 

Upper 

(MI)  Mitotic 
index  (%) 

Jejunum 

Savidge  et  al.,  1995 

0.97* 

1.05* 

(. PPG ) 

Proliferating  crypt 
cells  (cells/crypt) 

Jejunum 

Savidge  et  al.,  1995 

521.3* 

542.7* 

(. MPC )  Maturing 
crypt  cells 
(cells/crypt) 

Jejunum 

Savidge  et  al.,  1995 

114.7* 

119.3* 

(CPV)  Crypts  per 
villus 

Jejunum 

Ferguson  et  al.,  1977 

3.78*’+ 

6.5*’+ 

(' VH )  Villus  height 

(H 

Jejunum 

Smit  et  al.,  1986 

427.9* 

566.1* 

( VHC )  Villus  Duodenum 

height  to  cell  count 
conversion 
(//m/cells) 


(: Tm )  Length  of  M 
phase  (h) 

Duodenojejunum  Wright  et  al.,  1973 

1.09**  1.31** 

(Ts)  Length  of  S 
phase  (h) 

Jejunum  Schultze  et  al.,  1972 

O 

00 

o 

*Lower  or  upper 

bound  for  95%  CI. 

**Only  two  measurements  provided. 

1  Removed  an  outlier  of  15.4. 

°Measured  in  mice,  two  values  were  calculated  using  different  measuring  techniques. 


Hasan  and  Ferguson,  2.639*  3.679* 

1981 
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Table  3.4:  Human  parameter  associations  with  experimental  measurements. 


Parameter 

Biological  Meaning 

Formula  § 

7 

Steady-state  rate  of  proliferative  compartment  re¬ 
population 

(MI/Tm) 

a 

Maximum  rate  of  proliferative  compartment  re¬ 
population 

1  /(Tm  +  Ts) 

s 

Minimum  transit  rate  through  maturation  com¬ 

7*  PPC 

partment 

MPC(  1  + 1) 

Minimum  transit  rate  through  villus  compart¬ 

7  *  PPC  *  CPV 

ment 

VHC  *  VH(  1  +  to) 

§  Experimental  measurements:  mitosis  index  (MI),  proliferating  cells  per  crypt 
(. PPC ),  maturing  cells  per  crypt  ( MPC ),  crypts  per  villus  ( CPV ),  villus  height 
(’ VH ),  villus  height  to  cell  count  ( VHC ),  duration  of  M  phase  (Tm)  and  duration 
of  S  phase  (Tg). 
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3.5.2  Radiation  and  Burn  Data 

There  is  a  large  amount  of  time-dependent  acute  radiation  data  of  crypt  and  villus  cell  counts 
from  rodents  (Table  3.5).  We  assume  that  normalized  epithelial  cell  dynamics  of  mice  and 
rats  are  the  same  and  use  this  data  to  parameterize  our  murine  model.  In  each  study  that 
reports  the  radiation  dose  as  free-in-air  (FIA),  the  reported  doses  have  been  converted  to 
midline  tissue  (MLT)  in  Dant,  2016  using  radiation  and  animal  conversion  factors.  Our 
mathematical  model  simulates  MLT  radiation  doses  (Gy).  For  visualization  purposes,  we 
present  this  data  in  three  radiation  dose  ranges:  0  to  6  Gy  (Figure  3.3),  6  to  8  Gy  (Figure 
3.4),  and  8  to  25  Gy  (Figure  3.5). 


Table  3.5:  Optimization  data  for  the  irradiated  murine  small  intestine  model. 


Reference 

Measurements 

Species 

SI  Section 

MLT  Dose  (Gy) 

Galjaard  and  Bootsma,  1969 

C 

Rat 

Duodenum 

3.17 

Galjaard  et  al.,  1972 

C 

Rat 

Duodenum 

5.55 

Kononenko  and  Farafonov,  1969 

c,  V,  E 

Mouse 

ns 

7,  9 

Matsuzawa  and  Wilson,  1965 

C,  V 

Mouse 

Ileum 

26.14 

Matsuzawa  et  al.,  1973 

C,  V 

Mouse 

Duodenum 

4.34,  17.34 

Patel  et  al.,  2012 

C,  V 

Mouse 

Jejunum 

7.6 

Potten  and  Loeffler,  1990 

c 

Mouse 

Ileum 

6.08 

Quastler,  1956 

V 

Mouse 

ns 

1.74,  4.36,  17.43,  52.28 

Rijke,  1977 

C,  V 

Rat 

Duodenum 

2.38,  5.55 

Sato  et  al.,  1972 

C,  V 

Mouse 

Duodenum 

8.71 

Sherman  and  Quastler,  1960 

C,  V 

Mouse 

Ileum 

25 

Smirnova,  2010* 

C,  V 

Mouse 

ns 

7 

Note:  C-crypt  cells,  V-villus  cells,  E-total  epithelial  cells,  ns-not  specified. 
*  Original  reference  unknown. 


Crypt  cells  (0  to  6  Gy) 


"D 

d) 

N 

"nj 


1.2- 


o 

c 


1 0.8- 


o 

o 


Cl) 

CL 


^  0.4- 

a> 

O 


Experiment 

•  2.38  Gy  (Rijke  1977) 

A  3.17  Gy  (Galjaard  1969) 

■  4.34  Gy  (Matsuzawa  1 973) 
+  5.55  Gy  (Galjaard  1969) 

Bl  5.55  Gy  (Rijke  1977) 

2^5  5^0  7.5 

Time  after  irradiation  (days) 


Villus  cells  (0  to  6  Gy) 


"O 

CD 

N 

"o5 

E 

o 

c 


1.2 


fo. 

O 

o 

CD 

CL 

^0. 

0 

O 


• 

+ 

El 

•  ■ 

■ 

■ 

A 

• 

*  H 

■ 

* 

E 

E 

g|  M 

3 

Experiment 

•  1.74  Gy  (Quastlar  1956) 

A  2.38  Gy  (Rijke  1977) 

■  4.34  Gy  (Matsuzawa  1 973) 

+  4.36  Gy  (Quastlar  1956) 

El  5.55  Gy  (Rijke  1977) 

+ 

2.5  5.0  7.5 

Time  after  irradiation  (days) 


Figure  3.3:  Murine  data  from  acute  radiation  experiments  (0  to  6  Gy). 
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Figure  3.4:  Murine  data  from  acute  radiation  experiments  (6  to  8  Gy). 
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Figure  3.5:  Murine  data  from  acute  radiation  experiments  (8  to  25  Gy). 


Variability  in  data  must  be  considered  when  comparing  the  model  to  the  data.  Mea¬ 
surement  variability  can  be  expected  as  we  have  collected  data  from  studies  with  diverse 
experimental  procedures.  For  instance,  these  datasets  vary  in  histological  measurement 
techniques,  timing  of  measurements,  species  and  dose  estimation.  As  an  example  of  variabil¬ 
ity  in  histological  data,  Smirnova,  2010  reports  a  crypt  cell  nadir  approximately  70%  below 
normal  in  response  to  a  dose  of  7  Gy  (Figure  3.4),  whereas  no  more  than  a  50%  reduction  of 
crypt  cell  counts  are  reported  in  response  to  a  Dose  of  8.71  Gy  in  Sato  et  ah,  1972  (Figure 
3.5). 
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We  have  only  found  one  study  that  reports  time-dependent  human  epithelial  cell  counts 
after  radiation  exposure  (Trier  and  Browning,  1966).  In  Trier  and  Browning,  1966,  crypt 
epithelial  surface  layer  counts  are  reported  during,  and  following,  X-ray  therapy  for  a  clinical 
patient  with  seminoma  (a  germ  cell  tumor  of  the  testicle).  Unfortunately,  the  details  of  the 
timing  and  dosage  of  radiation  treatments  are  not  provided,  other  than  total  accumulated 
doses  approximately  1,  5,  8,  15,  22  and  28  days  after  the  start  of  treatment.  In  order  to 
recreate  the  experiment  in  silico,  we  estimate  the  dosing  regimen  by  assuming  equally  spaced 
treatments  between  each  reported  accumulated  dose  over  the  course  of  the  28-day  trial.  We 
have  also  converted  the  reported  doses  to  MLT  doses  using  a  0.96  conversion  factor  from 
roentgens  (R)  to  Gy,  and  a  0.667  conversion  factor  from  FIA  to  MLT  for  humans  (Anno 
et  ah,  2003).  Our  estimated  dosing  regimen  is  provided  in  Table  3.6. 


Table  3.6:  Estimated  dosing  regimen  for  the  Trier  and  Browning,  1966  experiment. 


Time  (d) 

FIA  Dose  (R) 

MLT  Dose  Received  (Gy) 

0 

100  (0  R  accumulated) 

0.64 

1.3 

200  (100  R  accumulated) 

1.28 

5.2 

116.67  (300  R  accumulated) 

0.747 

6.2 

116.67 

0.747 

7.1 

116.67 

0.747 

8.1 

143.75  (650  R  accumulated) 

0.92 

9.8 

143.75 

0.92 

10.7 

143.75 

0.92 

11.6 

143.75 

0.92 

12.5 

143.75 

0.92 

13.4 

143.75 

0.92 

14.2 

143.75 

0.92 

15.1 

133.33  (1800  R  accumulated) 

0.853 

16.3 

133.33 

0.853 

17.4 

133.33 

0.853 

18.6 

133.33 

0.853 

19.7 

133.33 

0.853 

20.8 

133.33 

0.853 

22 

140  (2600  R  accumulated) 

0.896 

23.2 

140 

0.896 

24.3 

140 

0.896 

25.5 

140 

0.896 

26.7 

140 

0.896 

28 

0  (3300  R  accumulated) 

0 

Lastly,  Table  3.7  lists  the  burn  and  combined  burn  and  radiation  experimental  data  we 
simulate  with  the  murine  model.  We  have  found  no  human  data  measuring  small  intestine 
epithelial  cellular  response  to  burn.  Due  to  a  limited  amount  of  data,  we  compare  the 
murine  model  to  normalized  mouse  and  rat  data,  taken  from  any  of  the  three  small  intestine 
segments  (ileum,  jejunum,  duodenum).  Although  the  model  was  developed  in  the  context  of 
the  murine  jejunum,  we  assume  the  normalized  kinetics  of  the  two  species  in  each  intestinal 
segment  are  the  same. 

A  limited  number  of  the  rodent  studies  reported  crypt  and  villus  cell  counts  following 
thermal  injury,  and  a  few  additional  studies  have  measured  villus  height  following  burn. 
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According  to  rat  data  from  Rijke,  1977,  normalized  villus  height  is  approximately  equal 
to  normalized  villus  cell  counts.  To  demonstrate  this,  we  present  data  from  Rijke,  1977 
in  Figure  B.13,  which  reported  villus  lengths  and  villus  cell  counts  after  various  radiation 
doses.  This  data  demonstrates  that  normalized  villus  height  is  a  reasonable  approximation 
of  normalized  villus  cell  counts. 

We  use  two  of  the  time-series  burn  datasets  to  fit  the  burn-related  parameters.  In  partic¬ 
ular,  we  use  the  villus  measurements  after  burn  from  Carter  et  ah,  2014  and  Jeschke  et  ah, 
2007.  Using  data  with  15%  and  60%  TBSA,  respectively  allows  the  model  to  develop  a 
dependence  on  %  TBSA.  The  remaining  burn,  radiation  and  combined  radiation  and  burn 
data  is  reserved  for  model  validation. 


Table  3.7:  Burn  data  simulated  by  the  murine  small  intestine  model. 


Reference 

Cells 

Data  Use 

Species 

SI  Section 

Insult 

Baker  and  Valeriote,  1968 

C 

Validation 

Rat 

D,  I 

16%  TBSA  +  7  Gy 

C 

Validation 

Rat 

D,  I 

7  Gy 

c 

Validation 

Rat 

D,  I 

16%  TBSA 

Carter  et  al.,  2014 

V 

Validation 

Mouse 

I 

15%  TBSA  +  4.72  Gy 

V 

Validation 

Mouse 

I 

4.72  Gy 

V 

Optimization 

Mouse 

I 

15%  TBSA 

Chen  et  al.,  2013 

V 

Validation 

Mouse 

I 

15%  TBSA 

Jeschke  et  al.,  2007 

V 

Optimization 

Rat 

All 

60%  TBSA 

Nesterenko  and  Piskarev,  1983* 

c 

Validation 

Rat 

na 

15%  TBSA  +  5.5  Gy 

c 

Validation 

Rat 

na 

5.5  Gy 

c 

Validation 

Rat 

na 

15%  TBSA 

Note:  C-crypt  cells,  V-villus  cells,  D-Duodenum,  I-Ileum,  na-not  available,  +  indicates  Cl. 
*Full  translation  unavailable. 
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3.5.3  Additional  Parameter  Bounds 


Using  the  methods  described  in  Section  2,  we  fit  all  parameters  of  the  murine  model  (exclud¬ 
ing  those  related  to  burn)  to  the  data  provided  in  Table  3.5.  When  fitting  the  parameters, 
we  restricted  7,  a,  5,  and  ^  to  the  narrow  regions  governed  by  the  values  and  equations 
presented  in  Tables  3.1  and  3.2,  respectively,  and  we  allowed  for  a  much  larger  range  of  mo¬ 
bility  for  the  remaining  parameters.  Table  3.8  lists  the  bounds  for  the  remaining  parameters 
enforced  when  exploring  parameter  space  for  improved  model  accuracy  of  both  human  and 
murine  models. 


Table  3.8:  Additional  parameter  bounds. 


Parameter 

Lower 

Upper 

9 

0.001 

1000 

0.001 

1000 

l 

1.5 

2.5 

m 

1.5 

2.5 

D0 

1  Gy 

10  Gy 

n 

1 

10 

V 

1  d"1 

10  d”1 

(f> 

0.001 

1000 

0.001 

1000 

CL\ 

0.001 

1000 

bo 

0.001 

1000 

bi 

0.001 

1000 

k 

0.001 

1000 

We  set  the  same  bounds  on  each  of  the  remaining  parameters  for  both  murine  and  human 
models.  We  restricted  6  and  t?  to  be  in  the  interval  [0.001,  1000],  which  allows  feedback  of 
maturing  crypt  and  villus  cells  to  be  up  to  three  orders  of  magnitude  less  than,  or  greater 
than,  feedback  provided  by  proliferating  cells.  I  and  m,  feedback  strengths  of  proliferation 
on  maturation  and  villus  transit  rates,  respectively,  are  each  set  to  2  in  Smirnova,  2010,  but 
we  do  not  have  access  to  the  referenced  study.  We  bound  l  and  m  in  the  interval  [1.5,  2.5], 
allowing  for  a  25%  adjustment  of  these  values. 

I11  Wentz,  2015,  the  multi-target,  single-hit  model  was  implemented  to  represent  dose 
response.  The  dose  response  parameters  (D0  and  n)  were  chosen  to  match  experimental 
data  approximated  from  Lesher  and  Lesher,  1970,  which  measured  the  number  of  crypt  cells 
in  synthesis  post-radiation.  Although  this  data  is  relevant,  irradiated  samples  are  rarely 
measured  at  the  same  time  as  controls,  which  makes  it  difficult  to  quantify  the  relative 
difference  between  experiment  and  control.  Instead  of  using  data  from  Lesher  and  Lesher, 
1970,  we  fit  the  dose  response  parameters  to  time-course  radiation  data  for  both  murine 
and  human  models,  while  placing  moderate  ranges  (within  one  order  of  magnitude  around 
the  values  chosen  in  Wentz,  2015)  on  both  parameters  (we  bound  Do  and  n  in  the  interval 
[1,10]). 
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The  remaining  radiation-related  parameters,  v  and  0,  determine  the  duration  and  strength 
by  which  the  damaged  cells  affect  proliferation,  respectively,  v  was  set  to  2d~1  chosen  in 
Smirnova,  2010,  assuming  average  apoptosis  lengths  of  12  hours.  In  Potten  and  Grant,  1998, 
it  was  determined  that  the  majority  of  apoptosis  occurs  within  3-6  hours  of  radiation,  sug¬ 
gesting  a  faster  rate  of  cell  death.  We  chose  to  bound  v  in  the  interval  [1, 10].  We  allow  0, 
the  relative  strength  of  feedback  of  the  damaged  cells  on  proliferation,  the  same  freedom  as 
9  and  d  ([0.001,1000]). 

The  burn-related  parameters,  do,  a±,  bo,  b\  and  k,  were  fit  to  the  rodent  villus  burn  data 
in  Carter  et  ah,  2014  and  Jeschke  et  ah,  2007.  Large  ranges  ([0.001,  1000])  were  allowed  for 
each  of  these  parameters  during  optimization,  as  no  prior  knowledge  was  available  for  these 
parameters.  Human  values  for  these  parameters  are  extrapolated  from  the  murine  model. 


3.6  Parameterization  Results 

We  now  present  model  simulations  of  the  experimental  data  used  to  optimize  the  radiation 
and  burn  kinetics  of  the  small  intestine  model.  Full  details  of  the  model  are  presented  in 
Appendix  A,  including  the  human  and  murine  parameters  (Table  A. 2).  Simulations  of  the 
models  are  presented  against  experimental  data  in  Figures  3.6  -  3.11.  When  available,  error 
bars  are  presented  as  ±  standard  deviations. 

In  addition  to  radiation  and  burn  insults,  we  also  simulated  rat  (Rijke  et  ah,  1976)  and 
mouse  (Wright  and  Al-Nafussi,  1982)  experiments,  presented  in  Figure  B.l.  In  Rijke  et 
ah,  1976,  a  method  called  villus  clamping  is  used  to  interrupt  the  blood  flow  to  the  small 
intestine,  causing  an  ischemic  villus  cell  damage,  killing  villus  cells  without  harming  crypt 
cells.  In  Wright  and  Al-Nafussi,  1982,  cytosine  arabinoside  is  used  to  kill  only  proliferating 
crypt  cells.  Both  of  these  experiments  measure  the  dynamics  of  small  intestinal  cell  recovery 
from  insults  other  than  radiation  or  burn. 

3.6.1  Murine  Radiation  Experiments 

Simulations  of  acute  radiation  rodent  experiments  are  presented  in  Figures  3.6  -  3.8.  In 
Section  3.5.2,  we  presented  this  data,  delineated  by  dose  and  cellular  measurement  (Figures 
3.3  -  3.5).  Considering  the  variability  in  the  data,  the  model  is  well  suited  for  capturing  each 
range  of  radiation  dose.  In  particular,  the  model  captures  the  time  and  value  of  villus  nadir 
in  a  dose-dependent  manner.  The  model  is  not  as  accurate  when  simulating  crypt  cells,  but 
it  is  possible  that  this  can  be  attributed  to  inconsistent  data. 
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Figure  3.6:  Simulations  of  a  murine  radiation  experiments  (0  to  6  Gy). 
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Figure  3.7:  Simulations  of  a  murine  radiation  experiments  (6  to  8  Gy). 
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Figure  3.8:  Simulations  of  a  murine  radiation  experiments  (8  to  25  Gy). 


22 


3.6.2  Human  Radiation  Experiment 

Trier  and  Browning,  1966  is  the  only  study  we  have  found  which  reports  time-dependent 
data  of  irradiated  human  epithelial  cells.  This  data  reports  villus  cell  counts  following  repeat 
exposures  of  an  unspecified  dosing  schedule,  which  we  estimated  (Table  3.6).  Figure  3.9 
provides  the  simulation  of  this  data  (red  line),  superimposed  on  the  data  (black  triangles). 
An  initial  drop  in  cell  counts  is  observed  in  the  data,  followed  by  a  quick  recovery  about 
12  days  after  the  beginning  of  the  treatment.  The  model  is  unable  to  capture  this  early 
recovery,  but  it  is  possible  that  this  phenomenon  is  an  artifact  of  measurement  error.  The 
model  accurately  captures  the  villus  nadir,  as  well  as  the  rate  of  recovery  following  the  end 
of  treatment. 

The  lack  of  adequate  human  data  is  the  biggest  limiting  factor  in  firmly  establishing 
human  parameters  for  the  model.  Currently  we  have  fit  7,  a,  6,  -0,  0,  d,  l ,  m,  Do,  n ,  v 
and  (f)  to  this  single  study,  although  7,  a,  5  and  0  have  been  restricted  to  small  ranges 
based  on  biological  measurements.  Allowing  this  many  parameters  to  depend  on  one  study 
of  one  individual  places  too  much  emphasis  on  this  data.  In  the  future,  if  we  are  unable  to 
find  more  data  of  this  nature,  we  will  explore  alternative  approaches  to  parameterizing  the 
human  model. 
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Figure  3.9:  Simulation  of  a  human  radiation  therapy  experiment. 
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3.6.3  Murine  Burn  Experiments 

In  this  section,  we  present  simulations  of  burn  response  in  the  rodent  small  intestine.  Figures 
3.10  and  3.11  provide  simulations  of  villus  and  crypt  cell  response  to  burn,  respectively.  In 
Chen  et  ah,  2013,  measurements  were  taken  after  two  different  drug  vehicles  were  applied 
(gavage  and  i.p.),  and  in  Baker  and  Valeriote,  1968,  measurements  were  taken  in  two  small 
intestine  segments.  The  burn  only  data  in  Carter  et  ah,  2014  and  Jeschke  et  ah,  2007  was 
used  to  optimize  the  burn-related  murine  parameters,  and  the  remaining  burn  data  serves 
to  validate  the  model. 
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Figure  3.10:  Simulations  of  villus  cells  from  rodent  burn  experiments. 
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Figure  3.11:  Simulations  of  crypt  cells  from  rodent  burn  experiments. 
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4  Combined  Injury 

In  the  event  of  a  nuclear  weapon  detonation,  it  is  anticipated  that  there  will  be  many 
casualties  exposed  to  various  ranges  of  acute  radiation,  burn,  trauma,  and  other  insults.  It 
is  therefore  important  that  our  models  capture  the  synergistic  effects  of  combined  injuries. 
In  this  section,  we  present  the  ability  of  the  mathematical  model  to  reproduce  combined 
radiation  and  burn  injury  experiments.  The  murine  model  is  validated  against  experimental 
data,  and  the  human  model  is  simulated  for  various  combinations  of  acute  radiation  and 
burn  injuries. 


4.1  Simulation  of  Murine  Data 

We  have  found  three  rodent  experiments  of  combined  radiation  and  burn  which  we  can  simu¬ 
late  with  our  model.  Each  of  these  experiments  measured  crypt  or  villus  cells  after  radiation, 
burn  and  combined  radiation  and  burn  injuries.  Figures  4.1  -  4.3  present  simulations  of  each 
of  these  experiments  against  the  experimental  data. 

In  Baker  and  Valeriote,  1968,  rat  duodenum  and  jejunum  crypt  cells  are  measured  four 
days  after  moderate  burn,  radiation  or  combined  insults  (Figure  4.1).  There  is  a  slight  de¬ 
crease  in  crypt  cells  four  days  after  radiation,  whereas  the  crypt  cell  population  has  recovered 
by  this  time  in  the  case  of  burn  and  combined  injuries.  The  model  is  able  to  capture  each 
of  these  features. 
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Figure  4.1:  Simulations  of  the  Baker  and  Valeriote,  1968  combined  injury  murine  experiment. 


The  Carter  et  ah,  2014  study  measured  villus  cells  1,  2  and  3  days  after  radiation,  burn 
and  combined  insults  (Figure  4.1).  The  model  captures  the  burn  response  of  this  experiment 
well,  which  is  to  be  expected  as  this  is  one  of  the  datasets  used  to  fit  the  model.  In  response 
to  radiation,  there  is  a  large  reduction  in  villus  cells  (more  than  30%),  the  extent  of  which 
is  not  reproduced  by  the  model.  Although  the  model  does  not  capture  the  quick  reduction 
in  villus  cells  reported  after  radiation,  the  model  is  able  to  reproduce  the  attenuation  of 
villus  cell  reduction  from  the  combined  injury.  In  particular,  radiation  and  burn  alone  result 
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in  an  approximate  20%  reduction  in  villus  cells  three  days  post-insult,  whereas  there  is  an 
approximate  40%  reduction  in  villus  cells  three  days  after  the  combined  injury. 


Figure  4.2:  Simulations  of  the  Carter  et  ah,  2014  combined  injury  murine  experiment. 


In  Nesterenko  and  Piskarev,  1983,  moderate  radiation,  burn  and  combined  insults  are 
applied  to  rats  before  crypt  cell  measurements  are  taken  2  hours  and  1,  3,  7,  14  and  30  days 
post-insult.  For  this  data,  we  have  not  been  able  to  acquire  a  full  translation  (Russian). 
It  is  possible  that  the  5.5  Gy  reported  in  this  experiment  is  FIA,  and  the  MLT  dose  is 
smaller.  In  this  case,  our  simulation  may  be  an  overestimate  of  the  response,  which  could 
account  for  the  large  villus  nadir  following  radiation  only.  The  granularity  in  the  reported 
measurements  could  also  account  for  the  differences  in  nadirs.  For  instance,  it  can  be  seen 
in  Figures  3. 3-3. 5  that  crypt  cells  reach  their  nadirs  approximately  two  days  after  radiation. 
This  could  account  for  differences  between  nadirs,  because  the  measurements  in  Nesterenko 
and  Piskarev,  1983  are  taken  at  1  and  3  days  post-radiation. 


Burn  (15%  TBSA)  Rad  (5.5  Gy)  Combined  Injury 


Figure  4.3:  Simulations  of  the  Nesterenko  and  Piskarev,  1983  combined  injury  murine  ex¬ 
periment. 
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4.2  Simulating  the  Human  Model 

We  have  extrapolated  human  values  of  the  burn-related  parameters  (do,  «i,  bo,  b\,  k  and  A&) 
from  the  murine  model.  Details  of  the  parameter  extrapolation  are  described  in  Section  C.  In 
short,  the  human  burn  parameters  were  fit  to  murine  burn  data,  scaled  in  time  and  amplitude 
based  on  quantitative  differences  in  human  and  murine  radiation  response.  In  Figures  4.4  and 
4.5,  we  provide  simulations  of  the  human  model  exposed  to  various  combinations  of  radiation 
and  burn  insults.  As  discussed  with  the  murine  model,  the  villus  nadir  has  very  little 
dependence  on  burn  size,  but  the  crypt  cell  proliferation  is  dampened.  This  human  model 
provides  simulations  of  human  response  to  combined  radiation  and  burn  injury.  Human 
burn  and  combined  injury  data  is  required  to  validate  this  model. 
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Figure  4.4:  Combined  injury  simulations  of  villus  cells  with  the  human  model. 
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Figure  4.5:  Combined  injury  simulations  of  crypt  cells  with  the  human  model. 
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5  Future  Work 

We  have  developed  mathematical  models  of  the  human  and  murine  small  intestine  epithelial 
cells  responding  to  radiation  and  burn  combined  injuries.  This  modeling  project  is  an  ongoing 
effort  that  can  be  improved  with  more  data  and  alternative  modeling  techniques.  This  section 
discusses  some  of  the  advancements  we  plan  to  make  to  this  model. 

5.1  Parameter  Analysis 

Parameter  sensitivity  analysis  and  uncertainty  analysis  are  important  tools  for  understanding 
the  inner  workings  of  a  model.  These  analysis  techniques  measure  the  dependencies  of  model 
outputs  to  parameter  errors  and  uncertainties.  A  handful  of  methods  have  been  developed  to 
accurately  and  efficiently  measure  parameter  sensitivities  and  uncertainties  (Hamby,  1994). 
In  the  future,  we  plan  to  leverage  these  methods  to  analyze  the  parameters  of  our  model. 
This  information  can  be  useful  for  detecting  model  error  and  examining  biological  anomalies. 

5.2  Alternative  Modeling  Approaches 

We  will  consider  modeling  other  animals  as  a  supplement  to  our  current  modeling  approach. 
For  instance,  non-human  primates  (NHPs)  have  been  studied  in  the  context  of  radiation- 
induced  gut  injuries  (MacVittie  et  ah,  2012b,  MacVittie  et  ah,  2012a,  Wang  et  ah,  2015). 
We  can  use  this  data  to  build  an  NHP  model  of  radiation  response  in  the  small  intestine. 
There  will  be  less  data  to  parameterize  an  NHP  model,  but  the  extrapolation  from  NHP  to 
human  would  be  more  reasonable  than  mice  to  humans. 

Citrulline  is  an  a-amino  acid  produced  mainly  by  enterocytes  in  the  small  intestine 
mucosa  and  has  been  identified  as  an  indicator  of  gut  damage  (Crenn  et  ah,  2000).  We  could 
extend  our  current  small  intestine  model  to  include  the  production  of  citrulline,  dependent 
on  the  number  of  epithelial  cells.  This  will  allow  us  to  leverage  currently  available  human 
citrulline  data.  Furthermore,  extending  the  small  intestine  model  to  include  citrulline  will 
increase  the  practicality  of  the  model.  Crypt  and  villus  measurements  require  invasive 
procedures,  whereas  citrulline  can  be  measured  through  plasma  extraction.  The  use  of 
citrulline  as  a  biomarker  can  be  beneficial  for  rapid  patient  analysis  for  injured  patients. 

5.3  Correlating  Model  Outputs  with  Clinical 
Endpoints 

An  upcoming  project  is  to  correlate  hemodynamic,  hematopoietic  and  small  intestine  kinetic 
models  to  clinical  endpoints.  For  instance,  we  plan  to  correlate  small  intestine  crypt  and 
villus  cell  counts,  following  radiation  and  burn  insults,  to  the  probability  of  mortality  and 
diarrhea.  Depending  on  the  available  data,  this  will  likely  include  measurements  of  crypt 
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and  villus  cell  nadirs  as  thresholds  for  probability  of  mortality.  There  may  be  more  high 
quality  citrulline  data,  in  which  case  we  may  reform  our  model  to  include  citrulline  produced 
from  epithelial  cells.  We  will  compare  our  results  to  other  models  that  have  already  made 
these  correlations  (Anno  et  ah,  1991,  for  example). 

5.4  Modeling  Additional  Insults 

In  addition  to  radiation  and  burn,  we  would  like  to  include  small  intestine  epithelial  cell 
response  to  various  types  of  trauma.  In  particular,  we  will  include  the  effects  of  hemorrhagic 
shock  on  the  intestinal  mucosa,  which  has  been  studied  in  rats  (Chang  et  ah,  2005).  Other 
effects  we  would  like  to  model  include  blast  overpressure,  penetrating  trauma  and  blunt 
trauma. 
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6  Conclusion 

The  small  intestine  is  essential  for  nutrient  absorption  and  is  vulnerable  to  insults  such 
as  thermal  injury  and  irradiation,  which  can  deteriorate  the  epithelial  lining  of  the  small 
intestine,  leading  to  increased  risk  of  internal  hemorrhaging,  bacterial  translocation,  sepsis 
and  mortality.  The  modeling  efforts  in  this  report  support  the  DTRA  mission,  which  aims  at 
providing  time-dependent  casualty  analysis  given  an  arbitrary  range,  and  combined  injuries 
expected  from  a  nuclear  attack,  ft  is  important  that  these  models  are  as  accurate  as  possible 
when  estimating  casualties  in  a  catastrophic  event. 

The  mathematical  models  developed  in  this  study  are  deterministic  sets  of  ODEs,  which 
capture  the  average  effects  murine  and  human  epithelial  cellular  response  to  radiation,  burn 
and  combined  injury.  The  murine  model  has  been  compared  to  mouse  and  rat  data,  collected 
from  each  of  the  three  sections  of  the  small  intestine.  Responding  to  radiation,  the  murine 
model  captures  the  average  of  crypt  and  villus  nadirs  (Figures  3.6  -  3.8),  which  provide 
biomarkers  for  epithelial  gut  damage.  We  will  also  explore  the  idea  of  using  a  time  element 
under  a  threshold  as  a  biomarker  for  intestinal  damage. 

The  human  model  is  limited  by  the  amount  of  data  available.  We  have  only  collected 
one  time-dependent  dataset  of  irradiated  human  epithelial  cells,  and  no  data  of  human 
epithelial  response  to  burn.  We  have  parameterized  the  effects  of  burn  in  the  human  model 
by  extrapolating  from  the  murine  model.  We  used  mouse  and  rat  data  to  parameterize 
the  murine  model,  and  the  murine  model  compares  well  to  available  burn,  and  combined 
radiation  and  burn  data.  The  aim  of  our  approach  is  to  improve  the  human  model  by 
supporting  it  with  a  well-developed  and  validated  murine  model.  We  hope  to  coordinate  with 
experimentalists  to  gather  more  accurate,  reliable  and  relevant  data  that  can  strengthen  our 
model.  Human  data  will  still  be  a  challenge  to  acquire  since  biopsies  of  the  small  intestine 
are  invasive  and  rarely  warranted  in  clinical  studies.  To  address  these  limitations,  we  will 
explore  alternative  biomarkers  of  intestinal  injury  such  as  plasma  citrullinc,  that  may  aid  in 
extrapolating  the  impact  of  radiation  on  the  GI  tract. 

In  this  study,  we  have  developed  a  mathematical  model  capable  of  simulating  combined 
radiation  and  burn  injury  response  of  small  intestine  cell  kinetics.  Data  restrictions  have 
provided  challenges  for  parameterizing  the  human  response  to  burn,  but  in  the  future  we 
aim  to  use  additional  approaches  for  parameter  calibration.  The  models  developed  in  this 
report  will  be  integrated  in  the  health  effects  modeling  software,  HENRE. 
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7  Abbreviations,  Acronyms  and  Symbols 


ARA 

Applied  Research  Associates,  Inc. 

ARS 

Acute  radiation  syndrome 

d 

Days 

DE 

Differential  Evolution 

DTRA 

Defense  Threat  Reduction  Agency 

FIA 

Free-in-Air 

FME 

Flexible  Modeling  Environment 

GIM 

Gut  Injury  Model 

Gy 

Gray 

HENRE 

Health  Effects  from  Nuclear  and  Radiological  Environments 

LPL 

Lethal  and  Potentially  Lethal 

MLT 

Midline  tissue 

na 

Not  available 

NHP 

Non-human  primate 

ODE 

Ordinary  differential  equation 

PAIR 

Proliferation  and  Intracellular  Repair 

RIPD 

Radiation-Induced  Performance  Decrement 

R 

Roentgen 

WMD 

Weapon  of  Mass  Destruction 

%  TBSA 

Percent  of  total  body  surface  area 
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A  Mathematical  Model  of  the  Small  Intestine 

Cell  Kinetics 

This  section  provides  a  comprehensive  list  of  parameters,  variables  and  equations  describing 
the  small  intestine  cell  kinetic  model.  Table  A.l  lists  each  parameter  and  variable  of  the 
model  with  their  biological  description. 


Table  A.l:  Biological  descriptions,  parameters  and  variables  in  the  small  intestine  mathe¬ 
matical  model. 


Component 

Biological  Description 

Units 

x  =  xud  +  xd 

Number  of  proliferating  crypt  cells 

cells/column 

j,ud 

Number  of  undamaged  proliferating  crypt  cells 

cells/column 

xd 

Number  of  damaged  proliferating  crypt  cells 

cells/column 

y 

Number  of  maturing  crypt  cells 

cells/column 

Z 

Number  of  villus  cells 

cells/column 

B 

Rate  of  proliferation 

d'1 

7 

Transit  rate  through  proliferative  compartment 

d"1 

F 

Transit  rate  through  maturation  compartment 

d"1 

E 

Transit  rate  through  villus  compartment 

d"1 

G 

Constant  upregulation  of  mediator 

(cells  /  column)  d“ 1 

H 

Rate  of  downregulation  of  mediator  by  epithelial  cells 

(cells  /  column)  “ 1  d~  1 

K 

Mediator  feedback  threshold 

cells/column 

a 

Maximum  rate  of  proliferative  compartment  repopulation 

d-1 

s 

Minimum  maturation  rate  through  maturation  compartment 

d-1 

Ip 

Minimum  transit  rate  through  villus  compartment 

d'1 

e 

Relative  effect  of  maturing  crypt  cells  on  repopulation  rate 

- 

i 9 

Relative  effect  of  villus  cells  on  repopulation  rate 

- 

L  =  lh 

Effect  of  proliferation  on  transit  rate  through  maturation  compartment 

d 

M  =  m/7 

Effect  of  proliferation  on  transit  rate  through  villus  compartment 

d 

D0 

Specifies  radiosensitivity  of  X  cells 

Gy 

n 

Number  of  targets,  X  cells 

- 

V 

Rate  of  damaged  cell  death 

d”1 

<t> 

Relative  effect  of  damaged  cells  on  repopulation  rate 

- 

s 

Burn  size 

%TBSA 

Go 

Determines  duration  of  burn  effect  on  villus  death 

d-1 

bo 

Describes  maximum  effect  on  villus  death 

- 

Cll 

Determines  duration  of  burn  effect  on  proliferation 

d-1 

hi 

Describes  maximum  effect  of  burn  on  proliferation 

- 

k 

Activation  threshold  for  burn  effect 

%TBSA 

A& 

Delay  after  burn  before  proliferation  suppression 

d 
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Equations  A.l  -  A. 8  define  the  normalized  small  intestine  cell  kinetic  model,  where  we 
denote  the  normalized  undamaged  proliferating  cells,  damaged  proliferating  cells,  maturing 
crypt  cells  and  villus  cells  as  xud,  xd,  y  and  z,  respectively.  Using  the  initial  conditions  in 
Equations  A. 9  and  A.  10,  the  model  simulates  an  acute  combined  radiation  dose  of  D  Gy 
(MLT)  and  thermal  injury  (S%  TBSA)  applied  at  time  t  =  0  d. 


A,ud 

X 

~d 

X 

y 


z 


B 


b 

W0(S,t ) 

WdS,t) 


7 xud(B  -  1)  -  W1(S,t)xud 

~d 

~VX 


5  ((1  +  l)xud  -  (1  +  lB)y) 

+  lB)y  -  (1  +  mB)ij  -  WQ(S,  t)z 
_ a/7 _ 

1  +  b(xud  +  (j)Xd  + 
a/7  —  1 


1  + 
b0 


O'y 

<5(1+0 

s 


+ 


#7 


S  +  k 


—aot 


b\  s+ke  ai^100  ^  t  >  Ab 
0  t<  Ab 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 
(A.  7) 
(A.8) 


x“d(0)  =  1  -  (1  -  e~D/D°)n  (A. 9) 

xd(0)  =  (l-e~D/D0)n  (A. 10) 

The  parameters  in  Table  A.2  can  be  used  to  simulate  small  intestine  cell  kinetics  of  either 
humans  or  mice. 
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Table  A. 2:  Murine  and  human  parameters. 


Parameter 

Murine  Value 

Human  Value 

7 

1.446  cl"1 

0.18  d”1 

a 

2.87  d"1 

2.589  d"1 

l 

2.456 

1.64 

8 

0.341  d"1 

0.305  d"1 

m 

1.639 

1.889 

ip 

0.284  d  1 

0.064  d"1 

e 

942.632 

11.899 

•& 

35.633 

4.633 

D0 

6.054  Gy 

6.005  Gy 

n 

1 

2 

V 

1.011  d"1 

2.727  d"1 

<P 

108.193 

743.469 

a0 

54.919  d”1 

1.458  d"1 

bo 

47.986 

1.346 

ax 

0.001  d"1 

0.001  d"1 

bx 

0.767 

0.123 

k 

0.008%  TBS  A 

1%  TBSA 

Ah 

1  d 

1  d 
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B  Supplementary  Figures 

This  section  provides  figures  supplemental  to  the  report.  Figure  B.l  presents  simulations  of 
epithelial  recovery  experiments  of  rodents  from  insults  other  than  radiation  or  burn.  Figures 
B.2-B.12  present  simulations  of  rodent  acute  radiation  experiments.  All  of  the  data  in  this 
section  were  used  to  optimize  parameters  of  the  murine  model.  In  each  figure,  simulation 
results  (solid  red  lines)  are  compared  to  experimental  data  (black  triangles),  and  error  bars 
represent  ±  standard  deviations. 

B.l  Additional  Rodent  Insult  Experiments 

We  examined  the  ability  of  the  model  to  simulate  crypt  and  villus  data  from  experiments 
in  Rijke  et  ah,  1976  and  Wright  and  Al-Nafussi,  1982.  In  Rijke  et  ah,  1976,  villus  clamping 
killed  villus  cells  in  the  small  intestine  of  a  rat,  and  crypt  and  villus  cells  were  measured 
2,  8,  16  and  24  hours  post-clamping.  In  Wright  and  Al-Nafussi,  1982,  cytosine  arabinoside 
killed  off  most  proliferating  crypt  cells  before  the  epithelial  cells  were  allowed  to  return  to 
their  original  values.  During  the  return  to  homeostasis,  crypt  and  villus  cells  were  measured 
once  daily  for  18  days  following  the  cytosine  arabinoside  injections.  Figure  B.l  shows  the 
model  simulations  of  these  experiments  with  the  murine  small  intestine  model.  Although 
the  model  initially  captures  the  dynamics  of  crypt  and  villus  cells  in  response  to  cytosine 
arabinoside,  the  amplitude  of  the  second  rebound  (negative  for  villus  cells,  positive  for  crypt 
cells)  is  not  reproduced. 
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Crypt  Cell  Mass  Crypt  Cell  Mass 

per  Column  (normalized)  per  Column  (normalized) 


Rijke  1976  -  Crypt  Cells 


Rijke  1976  -  Villus  Cells 


1.04- 


i.oo- 


0.96- 
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0  10  20 
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Wright  1982  -  Villus  Cells 


Time  after  injection  (days) 

Figure  B.l:  Simulations  of  non-radiation  rodent  insult  experiments. 
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B.2  Rodent  Radiation  Experiments 

The  following  section  provides  simulations  of  the  murine  model  against  mice  and  rat  acute 
radiation  experimental  data.  This  data  was  used  with  the  data  in  Section  B.l  to  optimize  all 
parameters  (not  related  to  burn)  of  the  murine  model.  The  complete  list  of  rodent  radiation 
data  can  be  found  in  Table  3.5. 


Galjaard  1969  Galjaard  1972 

Crypt  Cells  (3.1 7  Gy)  Crypt  Cells  (5.55  Gy) 


Figure  B.2:  Simulations  of  the  Galjaard  and  Bootsma,  1969  and  Galjaard  et  ah,  1972  ex¬ 
periments. 
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Figure  B.3:  Simulations  of  the  Kononenko  and  Farafonov,  1969  experiments. 
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Crypt  Cell  Mass 
per  Column  (normalized) 


Matsuzawa  1 965  Matsuzawa  1 965 

Crypt  Cells  (26.1 4  Gy)  Villus  Cells  (26.1 4  Gy) 


Figure  B.4:  Simulations  of  the  Matsuzawa  and  Wilson,  1965  experiments. 


45 


Crypt  Cell  Mass  Crypt  Cell  Mass 

per  Column  (normalized)  per  Column  (normalized) 
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Figure  B.5:  Simulations  of  the  Matsuzawa  et  al.,  1973  experiments. 
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Figure  B.6:  Simulations  of  the  Patel  et  ah,  2012  experiments. 


Potten  1 990 
Crypt  Cells  (6.08  Gy) 


Time  after  irradiation  (days) 

Figure  B.7:  Simulations  of  the  Potten  and  Loeffler,  1990  experiments. 
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Villus  Cells  Villus  Cells 

per  Column  (normalized)  per  Column  (normalized) 
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Figure  B.8:  Simulations  of  the  Quastler,  1956  experiments. 
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Crypt  Cell  Mass  Crypt  Cell  Mass 

per  Column  (normalized)  per  Column  (normalized) 
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Figure  B.9:  Simulations  of  the  Rijke,  1977  experiments. 
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Crypt  Cell  Mass  Crypt  Cell  Mass 

per  Column  (normalized)  per  Column  (normalized) 


Sato  1972 

Crypt  Cells  (8.71  Gy) 
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Figure  B.10:  Simulations  of  the  Sato  et  al.,  1972  experiments. 
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Figure  B.ll:  Simulations  of  the  Sherman  and  Quastlcr,  1960  experiment. 
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Crypt  Cell  Mass 
per  Column  (normalized) 
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Figure  B.12:  Simulations  of  an  experiment  reported  in  Smirnova,  2010. 
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B.3  Villus  Height  and  Cell  Count 

According  to  rat  data  from  Rijke,  1977,  normalized  villus  height  is  approximately  equal  to 
normalized  villus  cell  counts.  This  is  reasonable,  as  the  circumference  of  each  villus  remains 
relatively  constant  as  the  length  increases.  To  demonstrate  this,  we  present  data  from  Rijke, 
1977  in  Figure  B.13  (red  triangles),  which  reported  villus  lengths  and  villus  cell  counts  after 
various  radiation  doses.  The  best-fit  line  for  this  data  including  the  origin  (black  dashed 
line)  is  close  to  the  identity  (y  =  0.016  +  0.94  •  x)  and  fits  the  data  very  well  (r2  =  0.993). 
From  this,  we  assume  normalized  villus  height  is  a  reasonable  approximation  of  normalized 
villus  cell  counts. 
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Figure  B.13:  Plot  of  villus  height  against  villus  counts  measured  in  Rijke,  1977. 
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C  Human  Extrapolation 

We  do  not  have  the  human  data  required  to  quantify  the  components  of  burn  response,  which 
makes  it  very  difficult  to  extrapolate  parameters  for  a  human  model.  To  circumnavigate  this 
problem,  we  extrapolate  parameters  from  murine  kinetics.  For  our  approach,  we  compare 
human  and  murine  responses  to  radiation  in  order  to  make  assumptions  about  murine  and 
human  differences  in  burn  response. 

Two  key  components  of  response  to  insult  are  timing  and  amplitude.  We  fit  the  human 
model  to  villus  cell  data  which  has  been  scaled  (in  time  and  amplitude)  from  murine  burn 
data.  We  set  the  scaling  factors  based  on  the  timing  and  amplitude  of  human  and  murine 
responses  to  radiation.  In  particular,  we  compared  the  villus  cell  count  nadirs  following  an 
acute  radiation  dose  of  5  Gy.  With  these  relationships,  we  extrapolated  synthetic  human 
burn  data  from  the  data  used  to  optimize  the  murine  burn  parameters  (Table  3.7). 

In  this  section,  we  differentiate  between  murine  and  human  parameters  with  m  and  h 
subscripts,  respectively  (the  human  parameter  a 0  becomes  a0th,  for  example).  For  simplic¬ 
ity,  we  assume  the  delay  time  of  proliferation  suppression  is  equal  for  murine  and  human 
models  (A&  /,  =  Abim  =  1).  We  found  the  normalized  villus  cell  nadirs  of  the  murine  model 
(,tn5,m,  Zn5,m)  and  the  human  model  (tn5,/i,  zn§ ,h),  and  used  these  values  to  scale  the  time  and 
values  of  normalized  murine  villus  cell  burn  data  to  produce  synthetic  human  data: 


/  tn5,h 

'  f 

Ln5,m 


Vn5,h 
V nh,m 


Vi). 


(C.l) 


We  fit  the  human  burn  parameters  to  the  synthetic  data  (Eq.C.l),  adopting  the  same 
approach  used  to  parameterize  the  murine  burn  response.  The  parameter  values  resulting 
from  this  approach  are  provided  in  Table  C.l  and  simulations  of  moderate  radiation  (5  Gy) 
and  burn  (15%  TBSA)  insults  with  these  parameter  values  are  presented  in  Figures  C.l  and 
C.2. 


Table  C.l:  Extrapolated  human  burn  parameters 


Parameter 

Human  Value 

Murine  Value 

a0 

54.919  d’1 

1.458  d"1 

bo 

47.986 

1.346 

ai 

0.001  d"1 

0.001  d"1 

h 

0.767 

0.123 

k 

0.008%  TBSA 

1%  TBSA 

Ab 

1  d 

1  d 
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Human  Response  (Radiation) 


Time  after  radiation  (days) 


Human  Response  (Burn) 


Time  after  burn  (days) 


Murine  Response  (Radiation)  Murine  Response  (Burn) 


Figure  C.l:  Crypt  cell  simulations  for  human  and  murine  models  responding  to  radiation 
and  burn. 
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Human  Response  (Radiation) 


Time  after  radiation  (days) 


Human  Response  (Burn) 


Time  after  burn  (days) 


Murine  Response  (Radiation) 


Time  after  radiation  (days) 

Figure  C.2:  Villus  cell  simulations  for  human 
and  burn. 


Murine  Response  (Burn) 


Time  after  burn  (days) 

and  murine  models  responding  to  radiation 
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